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ABSTRACT 

Using an N -body evolution code that does not rely on softened potentials, I have created a suite of 
unbound interacting cluster pair simulations. The motions of the centers of mass of the clusters have 
been tracked and compared to the trajectories of point masses interacting via one of four different 
softened potential prescriptions. I find that the relationship between the impact parameter of the 
cluster interaction and the point-mass softening length that best approximates each cluster’s center- 
of-mass motion depends on the adopted prescription. In general, the range of allowed softening lengths 
grows roughly linearly with the impact parameter, but zero softening is acceptable in the majority 
of situations. In an IV-body simulation that adopts a fixed softening length, such relationships lead 
to the possibility of two-body effects, like dynamical friction, being either larger or smaller than 
the corresponding cluster situation. Further consideration of more specific IV-body situations leads 
estimating that a very small fraction of point-mass encounters experience two-body effects significantly 
different than those of equivalent clusters. 

Subject headings: methods:numerical — gravitation 


1. INTRODUCTION 

It is common for IV-body integration schemes to rely 
on softened potentials to speed up calculations. Small 
impact parameter, two-body encounters in a Newto¬ 
nian potential can drive timesteps to ve ry sma ll values 
in some IV-body implementations ( e.g lAarsethl 1200IL 
NBODY2). Softened potentials relax this behavior by 
capping the maximum acceleration any particle can ex¬ 
perience, thereby placing a lower limit on the timestep. 
From another perspective, softened potentials allow an 
IV-body system to behave more like a collisionless sys¬ 
tem, where two-body effects are absent. Minimizing 
strong scattering has led other IV-body codes to adopt 
softened potentials (e.a .. iSpringelH2005L GADGET-2). 

One way of bringing some physical reality to the idea 
of using softened potentials is to imagine that each par¬ 
ticle in a simulation is actually a cluster of unresolved 
particles. This work investigates the interactions of bi¬ 
nary clusters of particles relying strictly on Newtonian 
potentials (but note that throughout this work these 
clusters are not bound to one another). This is made 
possible by using the graphics processing uni t (GPU) 
enhan ced version of Aarseth’s NBODY6 code (lAarsethl 
120031 iNitadori fc Aarsethl I2012D . For comparison pur¬ 
poses, all simulations have been completed on a multi¬ 
core, 64-bit, dual-processor machine with an NVIDIA 
Quadro K2000 GPU. The local configuration used in this 
work utilizes GCC 4.4.7 and CUDA 4.0.8 compilers. 

I have focused on four different softened potential pre¬ 
scriptions to describe the two-body interactions that 
model the motions of each cluster’s center-of-mass. The 
first three belong to a class of prescription sometimes 
referred to as ‘Plummer-like’ softenings. These are the 
type of softening implemented in NBODY2. Prescription 
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one is described by 
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where r is the separation distance and f is the direction 
from one mass to the other. The gravitational constant 
G is set to unity in this work and the masses m\ = m 2 = 
1/2 when referring to the masses of clusters. Prescription 
two is described by 
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Note that this is not a useful potential, as it does not 
reduce to the Newtonian point-mass potential for 5 2 = 0. 
However, the force expression is Newtonian in that limit, 
allowing particle trajectories to be integrated. The third 
Plummer prescription is described by 
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The choice of these three prescriptions is not intended to 
be exhaustive but merely representative of the kinds of 
softenings one might adopt. 

The fourth prescription is a version of the ‘ spline ’ 
softening idea discussed by iHernauist fe Kal/j (|1989T ) 
in the context of smoothed particle hydrodynam- 
ics . Based on a smoothi ng ker nel function described 
in iMonaghan fe Lattanziol ( 19851 ). GADGET-2 smooths 
two-body interactions by assuming that the interacting 
particles are actually distributed masses “smeared” by 
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the smoothing kernel. The smoothing kernel compli¬ 
cates the potential and force expressions far beyond those 
above for the Plummer pres criptions; the details of their 
derivation may be found in iSpringel et all (l2001i l. Gen¬ 
erally, the spline softening approach uses three different 
approximations to the force, depending on the separa¬ 
tion of the particles. This behavior allows spline softened 
forces to exactly equal Newtonian values when particles 
are adequately separated. With Plummer softenings, the 
softened force is never exactly Newtonian, even though 
the difference becomes computationally negligible with 
adequate separation. In GADGET codes, the softening 
length is different than the smoothing parameter, and I 
have adopted the same multi plic ative f actor to connect 
them. As highlighted in ISpringel! (|2005T ) . this guarantees 
that the spline softened potential caused by a point mass 
when r = 0 is the same as that for a Plummer softened 
potential. 

Much of the basic physics of cluster-cluster interac¬ 
tions has been investigated in-depth over the past sev¬ 
eral decades. Early analytical attempts at describing the 
impact of a p assing mass o n a spherical system (e.g., 
iSnitzed 119581: lAlladinl 119651) inf ormed more n u merica l 
studies of the phenom ena (e.g., iRoos fc Normanl 119791 
lAguilar fe Whi tel 1 19851 and references therein). It is not 
the goal of this paper to revisit or duplicate the types 
of studies represented by these previous works. They are 
important to the present work only in the sense that they 
provide guideposts for the behavior of the clusters in the 
NBODY6 simulations. As each cluster’s center-of-mass in 
this work is intended to represent a point mass present in 
collisionless IV-body simulations, I focus on initial condi¬ 
tions that lead to hyperbolic orbits. The clusters here are 
intended to slide past one another without the possibility 
of becoming bound. I envision these clusters being what 
one would see if a microscope could be used to magnify 
the inner structure of individual masses in an A-body 
simulation that utilizes force softening. 

Previous discussions of softening in iV-body simula¬ 
tions can be grouped into different categories. Some re¬ 
searchers have focused on optimizing softening lengths 
to accurately reproduc e forces and/ or mi nimize comp u¬ 
tatio n al requiremen t s (lRqmeolll998t lAthanassoula et all 
120001 iPower et ali 120031) . Others have discussed 
more fundamental issues related to th e appropriateness 
of Plummer- like softening formulas (|Dver fc Id 119931 
IGerbeil 119961) or re-interpr eting so f tening prescriptions 
as smoothing operations dBarncsl 120121 ). The work 
presented here likewise investigates the impact that 
smoothed forces can have on the dynamics of an IV-body 
system. My addition to this chain is the investigation of 
how adequately softened point-mass motions can model 
extended cluster center-of-mass motions. In a practi¬ 
cal sense, the results of this work will hopefully inform 
the choices of future IV-body simulators. For example, 
if comparable dynamical accuracy can be achieved with 
smaller computational effort, perhaps a simple Plummer 
softening prescription would be more attractive than the 
(relatively) more complicated spline prescription. 

To begin, I present numerical details of the initializa¬ 
tion procedures (§ ED and the evolutions (§E2D- Dis¬ 
cussions of the main analyses of the dynamical evolutions 
and the key results that follow form the majority of Sec¬ 
tion [3] Finally, I present a summary of the techniques 


and outcomes in Section 2) 

2. NUMERICAL DETAILS 

2.1. Initial Conditions 

Clusters are composed of N c particles, where N c = 2 V 
and 10 < r] < 14, but rj = 13 (N c = 8192) has been 
adopted as the standard value. The positions of the par¬ 
ticles in one cluster are randomly chosen according to 
a spherical density distribution with maximum radius 
R c = 1. This work inve stigates clusters with equilibrium 
Plummer distributions (lPlummeiIll911f) with scalelength 
equal to one-tenth of R c and Gaussian distributions. The 
scalelength of the Gaussian distributions used is gener¬ 
ally set to R c /V 5; this makes the Gaussian clusters less 
centrally concentrated than the Plummer clusters. Par¬ 
ticles are given equal masses such that the total mass 
of a single cluster is M c = 1/2. For Plummer clusters, 
thermal velocities are isotropic and speeds are chosen 
from the Maxwellian speed distribution. An upper limit 
to particle speed is given by 3u r)rms at the appropriate 
radial location. With this choice, a single cluster is ini¬ 
tially closer to virial equilibrium than if t he upper limit 
were set to the escape speed (2'iy jrms , e.g.. lAarseth et aU 
11974 1. Particle velocities in Gaussian clusters are chosen 
to be isotropic and so that the cluster is in virial equi¬ 
librium. Dynamical evolutions of isolated, single clusters 
show a very modest mass loss; A M c /M c « 1 per cent for 
N c = 1024 and A M c /M c ss 0.5 per cent for N c = 16384. 
For single cluster NBODY6 evolutions, particles escape if 
they reach distances of 10 half-mass radii from the cen¬ 
ter of the cluster. The isolated cluster evolutions show 
evidence of some outer envelope expansion. The radius 
that contains 95 per cent of the mass of the cluster in¬ 
creases by approximately 5 per cent during an evolution 
lasting 10 single-cluster crossing times. For comparison, 
the radius containing 75 per cent of the mass increases 
by less than 1 per cent. 

Once a single cluster has been initialized, its center- 
of-mass location is moved from the origin and a mir¬ 
ror image of the cluster is created (different particle dis¬ 
tributions can be chosen for the two clusters; see Sec¬ 
tion [375]). Particle velocities are simply mirrored as well. 
Unless otherwise noted, initial center-of-mass separations 
in the x-direction Axhut are equal to four times the ini¬ 
tial cluster radius R Cl while there is no offset in the y- 
direction. The majority of simulations I have performed 
have center-of-mass separations in the z-direction (im¬ 
pact parameters) Azi n it that range between one and five 
times R c , but values up to 40i? c have been investigated. 

With all particle initial locations specified, the total 
potential energy 1/totai is determined. The self-potential 
energy of a cluster is the potential energy determined 
between all particles in that cluster. The total potential 
energy is the sum of the two cluster self-potential energies 
and an inter-cluster potential energy determined between 
particles in the different clusters. The speeds of each 
cluster’s center-of-mass are defined as fractions of J7 to tab 

^CM, initial = \J QUtotal, (4) 

where 0.2 < Q < 1.5 for the majority of the simula¬ 
tions discussed here. To provide some context to the 
Q values, I have compared the average thermal speed 
of particles in a Plummer cluster to the initial relative 
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center-of-mass speed for clusters. The average thermal 
speed is that of an equilibrium Plummer sphere and is in¬ 
variant to different cluster initial conditions. Varying the 
initial center-of-mass separation changes the inter-cluster 
potential energy, but by a very small amount relative to 
the self-potential energies. As a result, the total poten¬ 
tial energy is roughly equal to twice the self-potential of 
an equilibrium Plummer sphere, and the ratio of ther¬ 
mal to relative center-of-mass speeds is essentially just a 
function of Q, 

^CM, relative 
^thermal 

Specifically, Q values of 0.2, 0.5, 0.8, and 1.5 result in 
ratios of 1.0, 1.6, 2.0, and 2.7, respectively. The less- 
concentrated nature of the Gaussian clusters makes their 
potential energies smaller than those of Plummer clus¬ 
ters, with a corresponding reduction in the center-of- 
mass speeds of Gaussian versus Plummer clusters at a 
given Q. 

The initial center-of-mass velocities are chosen to lie in 
the x-direction. Finally, particles in each cluster have the 
appropriate center-of-mass velocity added to their pre¬ 
existing, thermal velocity. With all velocities assigned, 
an evolution timescale is calculated from, 


( 2M -) 5 ( 6 ) 

(2|F t otal|) 3 ’ ^ ' 

where F tota i is the sum of [/total, the center-of-mass ki¬ 
netic energies of the clusters, and the thermal kinetic 
energies of the clusters. 

2.2. Simulations 

The NBODY6 code advances all particles without soft¬ 
ening their interactions. All simulations have been 
evolved for a minimum of ten timescales. For Gaussian 
and Plummer clusters with Q = 0.2 (lower ucm, relative 
values), simulations have been extended to twenty and 
fifteen timescales, respectively. Visualizations of clus¬ 
ter interactions show that these values provide sufficient 
time for clusters to reach a point of minimum center- 
of-mass separation and then clearly separate. With the 
very different initial conditions involved in the binary 
cluster evolutions, the escaping particle criterion has 
been changed from the NBODY6 default behavior. For 
interacting cluster simulations, particles escape only if 
they reach distances 50 times the quantity (2M c ) 2 /[/ t otai. 
This allows simulations with large center-of-mass separa¬ 
tions (impact parameters up to 40 R c ) to be performed 
without artificially removing particles. With this new es¬ 
cape criterion, very few particles are lost from systems. 

During each evolution, particle positions and velocities 
are compiled every 0.1td yn . The total number of parti¬ 
cles in the system is reported at every output along with 
a list of particle identification numbers. By comparing 
lists of identifying particle numbers with the initial list, 
the cluster membership of each particle can be tracked. 
From this data, I calculate the centers-of-mass positions 
and velocities and then the thermal velocities of cluster 
particles. With those basics, I then calculate the ther¬ 
mal and center-of-mass kinetic energies of each cluster, 
the self-potential energies of each cluster, and the inter¬ 
cluster potential energy as functions of time. Total and 



2.2 y/Q. (5) 


individual cluster angular momenta and average cluster 
radius values are also tabulated as functions of time. 

The results of these simulations match those reported 
in earlier works deali ng wi t h similar system s and ini¬ 
tial conditions (e.a.. iDevadas Rao et a/.lll987lf . Interac¬ 
tions with impact parameters Azj n ; t > 2 R c are elastic 
in the sense that there is no change in either cluster’s 
center-of-mass kinetic energy during symmetric time in¬ 
tervals centered on the point of closest approach. For 
smaller impact parameters that lead to overlap of the 
systems at closest approach, the center-of-mass kinetic 
energies decrease. The self-potential energies of clus¬ 
ters generically increase. This agrees with the increase 
of average radius of each cluster during an interaction. 
Small impact parameter cases increase the average radii 
by roughly 15 per cent around the point of closest ap¬ 
proach, while larger impact parameters produce more 
modest changes, around 5 per cent. Since the clusters 
remain self-gravitating, their expansi on is coupled with 
a decrease in thermal kinetic energy dSpitzerl I1958T) . For 
larger impact parameter situations, the increase in self- 
potential has the same magnitude as the decrease in ther¬ 
mal kinetic energy. For smaller impact parameters, the 
loss in center-of-mass kinetic energy further boosts the 
self-potential change. 

Cluster angular momentum changes are generally very 
small. No angular momentum component changes by 
more than 1 per cent of the initial magnitude of the 
system angular momentum during the time interval sur¬ 
rounding the point of closest approach. Another per¬ 
spective on the size of these changes is that similarly 
sized changes occur when individual particles are lost 
from a cluster. Absent finite-particle-number variations, 
the clusters initially have zero angular momentum about 
their centers-of-mass. In the simulations with the small¬ 
est impact parameters, tidal interactions are able to 
torque them slightly. The small angular momenta ac¬ 
quired in these cases are oriented in the same direction 
as the initial total angular momentum. 

3. ANALYSIS & RESULTS 
3.1. Determining Softening Behaviors 

For any given set of initial conditions, the motions of 
cluster centers-of-mass have been matched to the mo¬ 
tions of point masses interacting via the softened poten¬ 
tials discussed in Section [T] Treating a softening length 
as a free pa rameter, the non -linear minimization rou¬ 
tine amoeba dPress et a/.lll994f) is used to determine what 
value of S produces the best representation of the center- 
of-mass motion. The figure of merit used by amoeba 
is a combination of least-squares deviations. The soft¬ 
ened motion starts from the simulated cluster initial con¬ 
ditions, and an adaptive timestep Runge-Kutta scheme 
(implemented in the IDL routine lsode) integrates soft¬ 
ened motion between the times at which simulated data 
(center-of-mass positions and velocities) exists. Given 
the M time values at which “data” values of x(ti), y(U), 
and z(ti) produced by NBODY6 exist, lsode advances 
the point mass locations and velocities from ti to U + i. 
The results shown and discussed here use a relative toler¬ 
ance of 10~ 7 ; fourth- and fifth-order Runge-Kutta steps 
result in relative differences smaller than the tolerance. 
Tests using a tolerance down to 10 -11 have also been per- 
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formed, but do not lead to any appreciable differences. 
Updating the time index i, the former final conditions 
become initial conditions and are recorded. In this way, 
a set of M corresponding softened “model” position and 
velocity values are produced. From these M pairs of val¬ 
ues, figures of merit that quantify differences between 
data and model values can be formed. The main figure 
of merit used in this work is formed by comparing posi¬ 
tion values. For each cluster, the M differences between 
center-of-mass and softened point x, y , and z positions 
are calculated. In each direction, and for each cluster, the 
average of the squared difference values is used to form 
an rms deviation. The two cluster rms values are then 
averaged in each direction. Finally, the figure of merit is 
formed by summing the directional rms values. In this 
way, the figure of merit reflects the average positional 
deviation of the model from the data over the entire sim¬ 
ulation. The optimal softening length 5 opt is that which 
produces the best agreement (smallest figure of merit) 
between the motion of the two clusters’ centers-of-mass 
and the softened point particles. 

With a dopt value, the figure-of-merit landscape is 
searched near the minimum to provide uncertainty val¬ 
ues. I have adopted a value of 10 -3 as the relevant figure- 
of-merit change, as that is roughly the size of the error in 
the location of a cluster center-of-mass. In sum, a given 
6 op t produces the best fit for a given model, and AS is 
determined by finding the two neighboring softening val¬ 
ues that produce figures of merit that are 10 -3 greater 
than the minimum value. Figure [1] shows a representa¬ 
tive example of the information used to determine un¬ 
certainties. In general, prescriptions one, two, and three 
form quadratic-like figure-of-merit curves. Prescription 
four tends to produce the kind of one-sided valleys seen 
in Figure [U This is unsurprising as the spline softening 
approach deals with ranges of softening values. What 
is mildly surprising is the fact that the prescription four 
£>o P t values always tend to be near the high-end of the 
uncertainty range. However, it is important to keep in 
mind that the difference in figure-of-merit values is ex¬ 
tremely small between <5 ~ 0 and (5 opt . In all but the 
smallest impact parameter situations, one can safely as¬ 
sume a softening length of zero with prescription four. 
For contrast with Figure [TJ Figure [2] focuses on the be¬ 
havior of the prescription four figure-of-merit curve for 
a simulation where Plummer clusters have Azi n it = Rc- 
Very small softening values are clearly ruled out in these 
kinds of interactions. 

This choice of figure of merit is straightforward but 
not unique. I have varied the calculation of <S op t to use 
differences between velocities rather than positions. Re¬ 
sulting values of <5 op t fall within the uncertainties derived 
from the position-based calculations. All remaining dis¬ 
cussions of <5 opt values and/or figure-of-merit values refer 
to the position-based approach. 

3.2. Comparison of Softening Prescriptions 

Two questions motivated this investigation. First, do 
particles interacting through softened potentials really 
behave like interacting clusters? Second, is any one pre¬ 
scription better than the others? I tackle the second 
question first in this section. However, note that this 
exercise is not determining an optimal form of softening 
potential, merely comparing several extant potentials. 


For any simulation, four <5 op t values are determined, 
one for each prescription. I use the corresponding best- 
fit figure-of-merit values as a comparison tool. For 25 
Plummer cluster simulations (N c = 8192, Aa"j rl j t = 4 R c , 
R c < Az; n it < 5 R c , 0.2 < Q < 5.0), prescription 1 pro¬ 
vides the best match 15 times, prescription 2 is best 4 
times, prescription 3 is best 1 time, and prescription 4 
is best 5 times. However, the differences between the 
values for the different prescriptions are generally very 
small. For example, in its 15 “wins”, the figure of merit 
for prescription 1 is, on average, smaller than its com¬ 
petitors by about 2 x 10 -4 . For the other prescriptions, 
their margins of victory are roughly 100 times smaller. 
Since these differences are smaller than the error in de¬ 
termining the center-of-mass position, I conclude that all 
of the prescriptions provide comparable levels of match. 

For a similar set of Gaussian cluster simulations, pre¬ 
scriptions 1 and 4 are evenly split; prescription 1 is best 
8 times, prescription 4 is best 10 times. The softness of 
the Gaussian clusters make it more difficult to determine 
how much better each prescription is. In simulations 
with Q < 0.5 and A 2 ; n ; t = i? c , the clusters experience 
strong tidal distortions. In the Q = 0.2 case, the two 
clusters actually merge during their encounter, making 
its analysis here worthless. Even though that case has 
been ignored here, it gives one the idea that figure-of- 
merit values will be generally larger than in the Plum¬ 
mer simulations, with correspondingly larger differences 
seen. In the 8 “wins” for prescription 1 (several of which 
are small impact parameter cases), the average figure 
of merit is about 0.01 smaller than that of its competi¬ 
tors. If the smallest impact parameter cases are removed, 
that value drops to about 10 -3 . Prescription 4 tends to 
win in higher impact parameter situations, so its margin 
of victory of 10 -5 is rather firm. While there is more 
hand-waving in this analysis, I again conclude that the 
prescriptions are basically equivalent in their ability to 
explain the center-of-mass motions of these clusters. 

3.3. Optimal Softening Lengths 

Figure [3] shows the results of determining the optimal 
softening length <5 op t using the four softening prescrip¬ 
tions given initially-equal-radius Plummer clusters with 
differing impact parameters in simulations where Q takes 
on the values 0.2, 0.5, 0.8, 1.5, and 5.0. A similar plot for 
Gaussian clusters is presented in Figure [4] Overall, the 
trend is for <5 op t values to decrease as the impact parame¬ 
ter increases. As discussed above, it is also typical to see 
non-zero lower limits to 5 for smaller impact parameters. 

In an attempt to filter out the impact of some other 
parameter values, I have created additional simulations 
with different initial separations and different cluster par¬ 
ticle numbers. Simulations with initial separations in 
the ^-direction twice the standard value produce S op t{b) 
curves comparable to those in Figure [3] Figure [5] shows 
how the S op t(b) curves differ between Plummer cluster 
Axi n it = 8 and Axi n i t = 4 simulations with N c = 8192. 
For clarity, only prescription one (5 opt values are shown. 
The level of agreement shown here is typical of what 
is seen for the other prescriptions and for results from 
Gaussian clusters. 

The results of varying particle numbers using Plum¬ 
mer clusters with Q = 0.5 are shown in Figure [Gj The 
panels isolate the different softening prescriptions, and 
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the different lines in each correspond to the various N c 
values. Overall, the clusters with the lowest resolution 
(N c = 1024) result in the largest values of (f op t- For 
N c > 2048, the values of <$ op t, uncertainties in d op t, and 
trends in d opt versus b are essentially the same. 

3.4. Results of Non-optimal Softening 

To investigate the impact of these S op t(b) relationships, 
imagine an TV-body simulation with a specified soften¬ 
ing length (5 s i m . Depending on the adopted softening 
prescription, one would like to choose <5 s ; m so that the 
majority of interactions would accurately describe the 
center-of-mass motion of two clusters of particles. 

For encounters with <5 s j m ^ S op t(b), the point masses 
will follow trajectories that do not match those of each 
cluster’s center-of-mass with the same impact parameter. 
To clarify the outcome of such encounters, further imag¬ 
ine simulating three pairs of point masses with the same 
impact parameter bo. In simulation A, <5 s ; m A > <5 0 pt(&o); 
in simulation B, <5 s ; m B < <5 0 pt(&o), while in simulation C, 
^simC = S op t(bo). I imagine the C pair as centers-of-mass 
for two clusters, since they follow the same trajectories. 
With a larger softening length, the pair in simulation A 
will experience a smaller change in velocity than those 
in simulation C. In general, point-mass interactions with 
Ssim > S op t(b) result in motions that underestimate the 
change in velocity relative to true cluster interactions. 
Conversely, the pair in simulation B will experience a 
larger change in velocity compared to those in simula¬ 
tion C; interactions with <f s i m < <5 opt (&) generically result 
in point-mass motions that overestimate the change in 
velocity relative to cluster interactions. If softened point- 
mass motions are intended to represent cluster center-of- 
mass motions, a single softening parameter is insufficient 
to correctly handle different impact parameters. 

It is standard to break the velocity change experienced 
by interacting point masses into components parallel and 
perpendicular to the initial direction of relative motion. 
The parallel component is related to the dynamical fric¬ 
tion a particle will experience while the perpendicular 
component relates to the relaxation time for the motion. 
In an iV-body simulation with a single fixed softening 
length 5 s i m , encounters where S s i m > (<)6 op t(&o) lead 
to weaker (stronger) dynamical friction and relaxation 
effects. Based on the ranges of (5 opt values presented 
above, very few situations will lead to an overestima¬ 
tion of dynamical friction/relaxation; S ss 0 is within the 
uncertainties in nearly all simulations. 

The question becomes, what is the likelihood for a par¬ 
ticle in a softened IV-body simulation to have an en¬ 
counter with any given impact parameter? I take the 
IV-body system to be a time-independent background 
volume density and investigate what a test particle mov¬ 
ing through that system would encounter. The IV-body 
system is spherical with radial extent TZ. For a test par¬ 
ticle on a radial orbit, I estimate the likelihood of an 
encounter by calculating the projected density of parti¬ 
cles in the system. This projected density is a function 
of a distance R measured in a plane perpendicular to the 
velocity of the particle on a radial orbit. The R value 
is just the impact parameter for the test particle and 
a ring of system particles centered on the origin. The 
projected density integrated around an annulus provides 
the number of system particles within a range of impact 


parameters. For a uniform density system, this annular 
density has a maximum value for R max /TZ =1/ V%- How¬ 
ever, more realistic centrally-concentrated systems have 
maximum annular densities that occur much closer to the 
center. As a simple example, consider a cored power-law 
volume density with asymptotic behaviors similar to the 
Plummer law, 

PO /yN 

P (1 + r/ro) 5 ’ [> 


where po is the maximum number density of system par¬ 
ticles and ro is a density profile scalelength. In this sit¬ 
uation, R max /TZ « 0.05 when TZ = 10r o . As a result, 
roughly 80 per cent of particles in such a system have 
b/TZ < 0.25 and about 15 per cent have b /TZ < 0.05. 
A Navarro-Frenk-White dNavarro et 7HfT997l NFW) vol¬ 
ume density produces a similar value of R max /TZ but has 
a more slowly decaying annulus density curve. Roughly 
80 per cent of particles in such a system have b/TZ < 0.55 
with about 10 per cent having b/TZ < 0.05. Steepening 
the outer volum e density profile to a Hernquist model 
(jHcrnquistl l l990 f shrinks R max /TZ to approximately 0.03 
and 80 per cent of particles have b/TZ < 0.35. Roughly 
10 per cent of particles have b/TZ < 0.03 in this case. 

For test particles on circular orbits, the approach is 
different. At any point along the orbit, the distribu¬ 
tion of impact parameters encountered by the particle 
is given by the projected background density function 
in the vicinity of the orbit. For a set of orbit radii 
0.01 < Ro < TZ, I have integrated the projected den¬ 
sity function over a circular region with radius b cen¬ 
tered on the orbit. For the most centrally-located cir¬ 
cular orbits investigated, the fraction of system particles 
with b < Rq grows to 10 per cent when Ro/TZ « 0.04 
for the Plummer-like cored power-law volume density. 
With a Hernquist system profile the values are very sim¬ 
ilar, but for an NFW volume density this level is not 
reached until Ro/TZ ss 0.08. For larger orbits with values 
of 0.1 < Ro/TZ < 0.3, roughly 25 per cent of system par¬ 
ticles have b < Rq, when the background is the Plummer- 
like density. This percentage increases to about 35 per 
cent for an NFW volume density. These sample calcula¬ 
tions suggest that in an isotropic IV-body system, with 
a mix of point masses on radial and circular orbits, it is 
not unreasonable to expect roughly 10 per cent of im¬ 
pact parameters to be small compared to the radius of 
the system. 

I now draw connections between the cluster radius R c , 
the size of a system in an iV-body simulation TZ, and 
the softening length. Taking the average density of one 
point mass spread over a cluster to be the same as the 
background system density p at a radius r, the cluster 
radius can be written as 


R 


3 

c 


3 

47 Tp(r) 


( 8 ) 


The minimum cluster radius is set by the maximum den¬ 
sity. Adopting the density of Equation [71 


R c 


3/ 

Nm 


1/3 

ro, 


(9) 


where A to t is the total number of particles in the simu- 
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lation and 



(1 + u ) 5 


drt. 


For 5 < TZ/ro < 50, I ft 0.08. This leads to 


distribution with a background density that is inversely 
proportional to radius produces no differential acceler¬ 
ations. Taking the background system to be composed 
of N other clusters, Ai = NM C . With the relationships 
developed earlier in this section, I estimate 


Rc « 0.17V- t 1/3 ^. (10) 

For an NFW background density, the right-hand side of 
Equation 1101 is approximately twice as big. As I am be¬ 
ing rather rough in this discussion, the factor of two is 
ignored and no distinction is made between the back¬ 
ground density profiles. 

For 10 5 < Wot < 10 7 , commonly realized values, 
Equation m leads to R r ft 1 0 3 1Z. Following the ar¬ 
gument in IPower et all (120031) . the IV-body simulation 
should not have softened accelerations larger than the 
mean-field acceleration at the system edge. This leads 
to the requirement that 


A > U 

^Sim _ rrrrz - 

(11) 

Combining Equations [10] and [TT| 


Rc < O.WsimA^t 6 . 

(12) 


For the same range of Wot mentioned above, EauationfT^l 
leads to 5 s \ m /R c > 1. 

Combined with the previous estimates of fractions of 
particles encountered by test particles, approximately 10 
per cent of interactions will have impact parameters less 
than about 40 R c . As this value is much larger than 
the range of impact parameters that have been discussed 
so far ( b/R c < 5), I have also run a limited number of 
binary cluster simulations with larger impact parameters 
(b/R c < 40). The 5 opt values for these simulations do not 
show monotonic trends, but the uncertainty ranges do. 
The lower limits of the ranges remain near zero, but the 
upper limits grow roughly linearly with b. 

For prescription one, the upper limits to 6 stay below 
R c until b ft 20i? c - For prescriptions two and three, 5/R c 
upper limits reach values of one at b ft 10i? c - With pre¬ 
scription four, the upper limit to S/R c ft 1 when b ft 2 R c . 
The inference I draw from these results is that adoption 
of the Power argument for choosing a simulation soften¬ 
ing length will lead to underestimating violent relaxation 
and relaxation effects in a very small percentage of inter¬ 
actions, those with very small impact parameters. 

One must be careful to not push this very simple anal¬ 
ysis too far. In an actual situation with clusters inter¬ 
acting as parts of a much larger system further compli¬ 
cations arise. In particular, tidal forces can affect the 
rates of separation of ea ch cluster’s ccntcr-of-mass (e.g., 
Privatika nto et al\ 120 1 61 Figure 1). A quick calculation 
suggests that the results presented here would not change 
dramatically in the presence of a background potential. 
Giving the background system a constant density (as in 
the core of an pseudo-isothermal halo) leads to a differ¬ 
ential gravitational acceleration field that is constant, 


d 9 


dr 


background 


GM 


(13) 


where A4 is the mass of the background system and 
G = 1 will be assumed from here on. I note that a cuspy 


_d£ 

dr 

where N = 10 6 has been adopted, in agreement with the 
other relationships. For comparison, I take the differen¬ 
tial gravitational acceleration due to a neighboring clus¬ 
ter to be that of a point mass (since the clusters should 
not overlap), 


M r 


background 


(low) 3 ’ 


(14) 


d g —2 M c 

dr cluster (aW) 3 ’ 


(15) 


where aR c is the distance from the center of the cluster. 
Taking a cluster to have these two contributions to tidal 
forces, we see that the relative strengths are comparable 
for a ft 10. As a < 10 at the points of closest approach 
for the simulations presented here, I estimate that back¬ 
ground tidal effects should be a small perturbation. Wit h 
an eye toward future work, note that iPriv atikant o et~cdi 
(|2016f) highlight that the details of the orientations of 
cluster orbits with respect to the larger system ( i.e ra¬ 
dially or tangentially biased) are important for accurate 
determination of tidal effects. Such an expansive inves¬ 
tigation is far beyond the scope of the present work. 


3.5. Repeated Interactions 

The expansion of clusters that accompanies interac¬ 
tions leads to the question of how cluster size affects 
the softening behavior. Specifically, if earlier interac¬ 
tions had expanded a cluster, will a subsequent interac¬ 
tion with an unexpanded cluster still be modeled by the 
same d op t(6) relationship? I note that adopting this as a 
physical picture requires one to think of clusters as self- 
gravitating systems that remain too small to be resolved 
in a larger-scale simulation. 

I have taken two approaches towards dealing with this 
situation. These simulations all start with the same 
center-of-mass initial conditions as the previously dis¬ 
cussed same-radius simulations. In the first approach, 
one of the clusters is allowed to have an outer radius 
twice as big as the other. The Plummer profiles of the 
two clusters are defined by the same scalelength, but 
one can be populated out to a larger radius. These 
different-radius-same-scalelength simulations model sit¬ 
uations where previous encounters have tidally extended 
the envelope of a cluster without affecting its core. The 
second approach assumes that the cluster core has also 
been strongly affected, with one cluster having twice the 
Plummer scalelength and outer radius of the other. 

Differences between the three families of simulations 
are evident in their energy and radius-change histories. 
Energies in the different-radius-same-scalelength simu¬ 
lations behave very much like those in the previously 
discussed same-radius cases (see Section 12.21) . In the 
different-radius-and-scalelength simulations, energy val¬ 
ues are now very different between the clusters, but the 
self-potential energy changes are still roughly of the same 
magnitude as the thermal kinetic energy changes in each 
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cluster. Clusters in all simulations expand as a result 
of their interactions. As with the same-radius situations, 
each cluster experiences a comparable fractional increase 
in radius in the different-radius-same-scalelength cases. 
For the different-radius-different-scalelengtli situations, 
the larger cluster experiences a smaller fractional in¬ 
crease in size compared to the small cluster but main¬ 
tains the largest size throughout the simulation. De¬ 
spite these differences, the 5 op t{b) relationships derived 
from the different-radius simulations are basically iden¬ 
tical to those from the same-radius simulations. Fig¬ 
ure [7] shows how the S op t{b) relationships for same- 
radius, different-radius-same-scalelength, and different- 
radius-different-scalelength simulations with N c = 8192 
compare. 

If the point masses of an ./V-body simulation were re¬ 
placed with clusters that have initially identical radii, 
encounters would quickly result in a distribution of clus¬ 
ter sizes and scalelengths. From the simple cases investi¬ 
gated here, it is reasonable to assume that the different 
cluster sizes will not appreciably alter the conclusions 
based on the same-size encounters. 

4. SUMMARY 

I have modeled the motions of interacting binary clus¬ 
ter centers-of-mass with a set of softened potential pre¬ 
scriptions. These kinds of simulations have a substan¬ 
tial history in the literature, where they have been used 
to investi gate galaxy-galaxy interactions in galaxy clus - 
ters (e.o.. lR.oos fc Normadl 197911 Aguilar fc Whitdil985l h 
The simulations discussed here are unique in that they 
cover a very different region of interaction parameter 
space and no softening is employed in the evolution of 
the clusters. 

The main inferences drawn from these simulations are: 

• Cluster centers-of-mass follow trajectories that can 
be well-described by point masses moving under 
the influence of softened potentials. 

• The softening prescriptions studied here provide 
matches to cluster center-of-mass motions that are 
comparable. Differences in the qualities of the fits 
provided by each prescription are generally smaller 
than the uncertainty in a cluster center-of-mass po¬ 
sition. 

• Optimal softening lengths can be found for any 
cluster-cluster interaction by minimizing discrep¬ 
ancies between cluster centers-of-mass trajectories 
and softened point-mass motions. For the Plum¬ 
mer prescriptions, optimal softening length values 
are largest for the closest encounters. For the spline 
prescription, the optimal softening length values 
grow nearly linearly with impact parameter. In¬ 
dependent of prescription, uncertainties in the op¬ 
timal values grow with the impact parameter. Soft¬ 
ening lengths near zero fall within the uncertainties 
for all but the closest encounters. 

• Given a prescription, the relationships between 
softening length and impact parameter for cluster 
interactions are robust against significant changes 
to initial conditions: cluster particle number, initial 


separation distance, relative velocity, and internal 
structure. 

• Softened point-mass interactions correctly describe 
cluster centers-of-mass motions when the softening 
length is optimal. For softenings less than the op¬ 
timal value, the point masses experience more re¬ 
laxation and dynamical friction effects than inter¬ 
acting clusters would. With softenings larger than 
optimal, point masses experience weaker two-body 
interactions than clusters would. 

• A negligible fraction of the interactions experienced 
by point masses in an ./V-body simulation will have 
impact parameters that place them in the regime of 
experiencing weaker relaxation and dynamical fric¬ 
tion effects than would clusters. However, collapse 
simulations that involve varying numbers of parti¬ 
cles with given impact parameters will most likely 
be affected differently, possibly more severely, than 
the static situations discussed here. 
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Figure 1. Panels a, b, c, and d show figure-of-merit curves for softening prescriptions one, two, three, and four, respectively. These curves 
are based on a Plummer cluster simulation with N c = 8192, Q = 0.2, Axi n it = 4 R c , Azinit = 2 R c . The thin horizontal line indicates the 
limiting figure-of-merit value, minimum plus 10“ 3 . The thick, solid vertical lines mark the low and high range values for the softening 
value, while the thick, dashed vertical line indicates the 5 op t value determined by the amoeba routine. Due to the rather large range 
imposed by the prescription 4 5 op t value and the rather small <5 op t value for prescription 1, it is nearly impossible to see the prescription 
1 curve. However, viewing that curve individually reveals it to behave very similarly to the prescription 2 and 3 curves. Note that for 
prescription 4 (panel d) the minimum is actually nearer the high end of the range. This is a generic behavior for prescription 4 (see, for 
example, Figure [31. 
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<5/R c 


Figure 2. A single, prescription four figure-of-merit curve for a Plummer cluster interaction with N c = 8192, Q = 0.5, Arrinit = 4 R c , and 
A^init = Rc- Contrast the quadratic-like behavior here with the one-sided trough seen in Figure [Til. Interactions with the smallest impact 
parameters investigated here tend to enforce lower limits to the softening length regardless of other initial conditions or cluster density 
profile. 
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Figure 3. The relationships between the optimal softening length 5 op t and impact parameter for the four softening prescriptions in 
simulations of Plummer clusters with N c = 8192. Lines connect the <5 op t values, while the errorbars indicate the ranges for those values. 
Each line represents results from simulations with different Q values. Panels a, b, c, and d correspond to prescriptions one, two, three, 
and four, respectively. Softening lengths and impact parameters are scaled by the initial cluster radius R c . Uncertainty ranges grow with 
impact parameter. For all but the smallest impact parameter cases, softening lengths of zero are reasonable. 
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12 3 4 5 

b/R c 

Figure 4. The relationships between the optimal softening length 6 op t and impact parameter for the four softening prescriptions in 
simulations of Gaussian clusters with N c = 8192. Panels and linestyles are the same as those in Figure l3l The Q = 0.2, b/R c = 1 results 
are not reliable as those simulations lead to cluster convergence. There is a trend for smaller impact parameter cases to result in larger 
5 op t values. As with the Plummer clusters, softening lengths of zero are reasonable for larger impact parameter cases. 
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Figure 5. Comparisons between prescription one softening lengths given different initial separations in the ^-direction. Panels a, b, c, 
d, and e correspond to results from simulations with Q = 0.2, Q = 0.5, Q = 0.8, Q = 1.5, and Q = 5.0, respectively. All simulations 
have N c = 8192. Versions of this plot using information from different prescriptions and/or cluster density profile show similar levels of 
agreement. 
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Figure 6. Comparisons between softening values given different cluster particle numbers. Panels a, b, c, and d correspond to the different 
softening prescriptions. The various linestyles correspond to different particle numbers 1024 < N c < 16384. In general, <5 op t values and 
uncertainties are consistent as long as N c > 2048. 
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Figure 7. Comparisons between S 0 pt(b) relationships for N c = 8192, Q = 0.5, A:ri n i t = 4 R c simulations in which the Plummer clusters 
have initially (i) same outer radii (‘same-size’ line), (ii) different outer radii but same density scalelengths (‘diff-size-same-scale’ line), and 
(iii) different outer radii and different Plummer density scalelengths (‘diff-size-diff-scale’ line). No different-size simulations have been 
created with b = 1, as they lead to excessive cluster overlap. In simulations with differences, the length in one cluster is twice that in the 
other cluster. For example, in the ‘diff-size-same-scale’ simulations, one cluster has twice the outer radius of the other. Panels a, b, c, 
and d correspond to prescriptions one, two, three, and four, respectively. In general, interactions with a “puffed up” cluster (larger outer 
radius and larger scalelength) produce larger <5 op t values. However, differences in 5 op t values are smaller than the uncertainties in all but 
the smallest impact parameter cases. 





































